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Abstract: We present an analytic technique for evaluating single cuts for one-loop in- 
tegrands, where exactly one propagator is taken to be on shell. Our method extends the 
double-cut integration formalism of one-loop amplitudes to the single-cut case. We argue 
that single cuts give meaningful information about amplitudes when taken at the integrand 
level. We discuss applications to the computation of tadpole coefficients. 
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1. Introduction 

Precision calculations in gauge theories such as QCD are needed for observations in hadron 
collider experiments and are an important motivation for studies of higher-order scattering 
in general. Recently there have been rapid developments in computational techniques, 
many centered around so-called unitarity methods [1,2]. 

The idea underlying the unitarity methods is to constrain amplitudes by their branch 
cuts in various channels. The constraints from cuts are sufficient only for certain classes 
of amplitudes, such as one-loop amplitudes in massless theories [3]. Unitarity cuts are 
evaluated by the Cutkosky rules [4], which put two propagators on shell. Rather than 
completing the dispersion integral at this point, unitarity methods store the imaginary 
part of the amplitude obtained this way, and accumulate similar information from all 
possible cuts in analytic continuations to other kinematic regions. 

Putting more than two propagators of a one-loop amplitude on shell gives a "general- 
ized" unitarity cut [5-16]. It isolates the part of the amplitude from diagrams containing 
the specific propagators being cut. An advantage of such cuts is that the delta functions 
that are used to place the propagators on shell effectively reduce the dimension of the 
remaining integral. So these generalized cuts are relatively easy to evaluate but give infor- 
mation about a correspondingly smaller part of the full amplitude. However, generalized 
unitarity can also be taken in the other direction, by cutting just one propagator. This is 
the "single cut" we discuss in this paper. Our focus differs from previous studies [17-20] 
in that we are interested in the explicit evaluation of the 4-dimensional single-cut integral, 
with standard Feynman propagators, with a view towards computing tadpole coefficients 
in amplitudes with internal masses. 
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One-loop amplitudes have an expansion in terms of master integrals, which are scalar 
box, triangle, bubble and tadpole integrals, with coefficients that are rational functions 
of the kinematic invariants [21-36]. In practice, unitarity methods operate by pattern 
matching cuts of amplitudes with cuts of master integrals. Tadpole integrals lack cuts in 
physical channels (double cuts), but we hope that single cuts should be useful to calculate 
their coefficients. In principle, single cuts should give information about all parts of the 
amplitude. We study the single cuts of master integrals and find that the tadpole gives a 
rational value, while the others have purely logarithmic single cuts. Therefore the tadpole 
coefficients can be targeted by discarding all logarithms. Other approaches to analytic com- 
putation of tadpole coefficients have been proposed, using universal UV and IR divergent 
behavior [37, 38] or introducing an auxiliary propagator [39] . 

Inspired by the formalism for explicit evaluation of double cuts [40-48], where the 
two-dimensional (or (D — 2)-dimensional) integral is performed algebraically by the Cauchy 
residue theorem, we rewrite our loop momentum in terms of spinor variables, which can in 
turn be exchanged for a complex variable and a real parameter. Since we now have fewer 
cut constraints than in a double cut, we are not able to evaluate the integral over the real 
parameter trivially. In fact, we leave it unevaluated. We find that it is sufficient to work 
at the integrand level; moreover, the full single cut integral will typically diverge. (Su- 
persymmetric theories are a notable exception: there, single cuts are well defined [18] and 
can be used to check expressions for planar multi-loop amplitudes [49].) The integral over 
the complex variable is addressed by the Generalized Cauchy Theorem, invoked similarly 
in [48]. If we are computing single cuts to find tadpole coefficients, we see another major 
difference compared to double cuts. While double cut evaluations were a matter of evalu- 
ating residues at poles, here it is the contour integral part of the formula that dominates 
the tadpole contribution, so we do not compute residues at all. 

Working at the integrand level is important because functions of the loop momentum 
can have non- vanishing single cuts even if they integrate to zero [18]. Therefore, a proper 
expansion of the integrand from a Feynman diagram includes all such terms, which have 
been thoroughly classified in the context of four-dimensional reduction [10,50]. 

The paper is organized as follows. In Section 2, we explain the phase space integration 
for single cuts and evaluate the cuts for master integrals. In Section 3, we present the single 
cuts of the first few types of integrands with tensor numerators. In Section 4, we extract 
the tadpole coefficients in some sample tensor integrands. In Section 5, we discuss the 
modifications necessary in the case of massless external legs (which are vanishing Gram 
determinants for bubble integrals). The appendix contains further details on the phase 
space parametrization. 

2. Evaluation of the single cut 

Our starting point is the one-loop integrand, 

j _ N(k) 

D D 1 ■ ■ ■ D k ' 
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(2.1) 



where N(k) is a polynomial in the loop momentum k, and the denominator factors are 



Di = (k-KiY-mf. (2.2) 

The single cut is a singularity of the amplitude selecting single propagators. We define the 
4-dimensional single-cut operator for a particular propagator Di to act on the integrand as 

A^/A^tA) (__^__). (,3) 

The single cut must be applied to the integrand, because there are non-vanishing contri- 
butions from so-called spurious terms (terms that vanish upon integration). Working with 
the integrand allows us to identify the particular propagator being cut. 

To evaluate the single cut analytically we introduce convenient variables and reference 
vectors following [42,48,51]. First, we would like to exchange the original loop momentum 
variable k for a null vector i\ in order to make use of the spinor formalism [52-56]. With 
respect to an arbitrary Lorentz vector K satisfying K 2 ^ and Kq > 0, we decompose the 
loop momentum variable, 

k = £i + £K, (2.4) 

where £ is a scalar factor taking a value such that i\ is null. This condition is implemented 
by another delta function and an integral over £, 



J d 4 k (.) = J di J d 4 h 5 (+) (^)(2£i • K)(t 



Let us now look at the single cut integral, Q2.3[) . For simplicity, suppose we have redefined 
the loop momentum so that Ki = 0. We comment further on allowable redefinitions below. 
Also, let m = 771$. Applying the change of variables ([2.4|), we have 

d'k 5^ (k 2 - m 2 )(.) = Jdi J d*h 5( + \£l)?^5(£ -£)(.), (2.5) 

where 

a = (« 1 . J o'+«w, i = - 2t \ K K t^ - 

We continue by exchanging the integration of i\ over the lightcone for an integral over the 
complex plane and a real parameter t. First, we express the vector K as a sum of two 
null vectors, p and q, which will later allow us to integrate over familiar complex variables 
instead of spinors. Since our choice of K was arbitrary, we can consider the choice of null 
p and q to be the starting point, making sure that p • q / and Po + qo > 0. 

K=p + q, p 2 = q 2 = 0. 
The replacement for the loop momentum is 

< =t(f + zzq» + I (q\^\p] - I <p|V*|g]) . (2.6) 
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Since 2p ■ q = K and 1l\ ■ K = t(l + zz)K , the integral measure becomes 



/ 



-*,»)(.)= ids r d 4 f{ — ;Jwp+^*K-e> (2 . 7) 



where now 



and we have defined 



2 ^/t 2 (l + zz) 2 + u + i(l + zz) 



(21 



4m 

u \ 



K 2 



We will find it convenient to work in the limit u — > 0, equivalent to choosing our arbitrary 
K such that K 2 3> m 2 . In fact, we will be able to set u = exactly in all the cuts we 
study. 

Now, the £ integral can be performed immediately by the delta function substitution. 
The complex integration over z and z will be performed by the Generalized Cauchy Formula 
as described in [48]. That is, for the integrand F(z,z) we construct a primitive G(z, z) 
with respect to, say, z. Let D be a disk in the complex plane encompassing all poles of 
G(z, z) viewed as a function of z. Then 

/ F(z,z) dz Adz = <p dz G(z,z) -2-iri V Res{G(z, z), Zj ). (2.9) 

JD JdD , 

poles Zj 

In practice, we use A to denote the radius of the disk D and rewrite the complex 
variable in terms of polar coordinates: 

z = re ia - D = {(r, a) \ < r < A, < a < 2vr}. 

Additional details on the phase space parametrization are given in the appendix. 

The final integration over t will not actually be carried out; we extract the information 
we need at the integrand level. In fact, the integral over t typically diverges. It should be 
cut off both from above and below, and any further transformations should be consistent 
with such a cutoff. For example, we allow linear shifts of the original loop momentum, but 
not global rescalings. 

Analytically, the t-dependence of the cut integrands turns out to be useful, as we 
can restrict our attention to leading or subleading terms in t in order to compute tadpole 
coefficients. 

In our first change of variables for the single cut, ( |2.5| ), we assumed that the cut 
propagator was already in the form (k 2 — m 2 ). In fact, it is good to redefine the loop 
momentum k — > k + Ki in each term of the single cut operation (^^) so that the delta 
function is 5^ + \k 2 — m 2 ). Because the single cut integral diverges, this redefinition is most 
obviously valid in the limit of large K 2 , which we will implement routinely. Otherwise, we 
must pay attention to the details of the cutoffs. 

We will now study the single cuts of the integrands of master integrals. 
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Single cut of tadpole 

The scalar tadpole is the simplest integrand allowing the single cut operation. With the 
variables of (|2.7p , 

dt f ,., , N K 2 t 2 (l + zz) 



ADo k 2 -m 2 
One z-primitive of the integrand is 



— / (idz A dz) 



v / t 2 (l + zz) 2 + 



Ti = {l/z)K 2 v / t 2 {l + zz) 2 + u, (2.10) 
so the result of applying the Generalized Cauchy Formula (|2.9|) is 



(2tt)K 2 (y/t 2 (l + A 2 ) 2 + u- V^T 



u 

The factor of i has dropped out in converting the differential form (see the appendix for 
details). 

In the limit of vanishing u, the single cut of the tadpole is therefore 

(2ir)K 2 t (a 2 + „, - <9N - A;) ■ 
v ' V 2t 2 (l + A 2 ) 2t 2 J 

We omit writing the integral over t (and the factor of 4 in the denominator), since it will 
not be necessary in this paper to carry out this integration. In subsequent integrals we will 
also drop all subleading terms in it so that the formulas are more manageable, while still 
sufficiently distinct. 

Single cut of bubble 

Consider the integrand 1/{DqD{), and take the single cut of Dq. 

1 



A 



D 



d ^ k 6 (+)(k 2 -m 2 )^- 



(k 2 - m 2 ){(k - K^ 2 - m 2 ) 
Now, under the integral with the delta function, 

D l = f l -2i l -K l -2iK-K u (2.11) 

where 

h = K 2 -m 2 + m 2 . (2.12) 

In the limit of vanishing u, we also have £ — > 0, so we can neglect the last term of ( |2. 1 1[ ) . 1 

, followed by expansion in the 



So we proceed with the replacement D\ = fx — t (^X\Ki\\ 
null vectors p and q. The single cut is 

K 2 t 

(idz Adz)—— -, (2.13) 

J-{Ki,z, z) 

where we have defined 

T(Ki,x, z) = h-t {{jp\Ki\p\ +z(p\K i \q]+z {q\Ki\p] + zz (q\Ki\q]) . (2.14) 
One z-primitive of the integrand in ( |2,13j ) is given by 

= K 2 log TjKuz^z) 
' 2 ~ (p\Ki\q\ +z(q\K 1 \q]' 1 ' > 



1 Indeed, one can check using (|2.q) that in the rest frame K = (Kq, 0, 0, 0), this term falls off as 1/Kq. 
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Single cut of triangle and box 

The analysis of triangle and box integrands is similar to the bubble. For the triangle 
integrand, 1/(DqD\D 2 ), the single cut can be expressed in terms of the primitive, given by 

V3= I d ~ Z HKuz^HK^z) (2 ' 16) 
^ lo g (5^4Y (2-17) 



V(K u K 2 ,z) *\F(K 2 ,z,z 

For the box 1 / \DqDiD 2 D 3 ) , the z-primitive is given by 

f KH 
J F{Ki,z,z)F{K 2 ,z,z)F{K 3 ,z,z) 
_ K 2 {\\Kx\q]\ogF{Ki,z,z) K 2 {X\K 2 \q} log F(K 2 , z, z) 



*D{Ki,K 2 , z)V(K 3 , Ki, z) V{K 2 ,K 3 ,z)V{K u K 2 ,z) 

K 2 (X\K 3 \q]logF(K 3 ,z,z) 
^ V{K 3 ,K x ,z)V{K 2 ,K 3 ,z) ' K - ' 

We have defined 

V{K h Kj,z) = fj (\\Ki\q] - fi (\\Kj\q] - t \pq] {X\KiKj\X) , where \\) = \p) + z \q) . 
Overview: master integrands and strategy for single cuts 

We have found the primitives V\,V2,V 3 ,V±, which are associated to the various master 
integrands. The single-cut calculation would be completed by the Generalized Cauchy 



formula, (2.9). It turns out that in the limit A — > 00, we can ignore the residues and 



restrict our attention to the closed line integral term. In polar coordinates, 

<b V n (z,z)dz = i da Ae ia V n (Ae ia ,Ae- ia ). (2.20) 

JdD JO 

In the limit A — > 00, the leading behavior of the integrands is 

Ae ia V! ~ A 2 t, (2.21) 



logfA 2 ) 

Ae ^n-^jT> n = 2,3,4. (2.22) 

The higher-point integrands are suppressed by powers of A. Moreover, among all these 
primitives, all rational terms come from tadpoles. Tadpole primitives are purely rational, 
while the others are purely logarithmic. Therefore, in an algorithm targeting tadpole 
coefficients, we will select terms of the single cut with specific dependence on A 2 . 

For fully generic integrands, we only need the single-cut operator selecting the terms 
proportional to A 2 t, which we denote by a bar. 

A Di [/] == A Di [/] . (2.23) 

A^t— terms 
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For integrands with null external momenta, it is convenient to define some further refine- 
ments of the operator, as follows. 



A Do [I] = A Do [I] 
A Do [I] = A Da [I] 



A 4 f 2 — terms 
\og((q\Ki \q] A 2 )— terms 



(2.24) 



In the case where the masses of the propagators are not all distinct, we may also want to 
collect the single cuts of all propagators with a given mass, as indicated by the subscript 



m 



A m2 [/] = A A W • ( 2 - 25 ) 

i:m,i=m 

It is now clear how to distinguish single cuts of the various master integrands. Since 
spurious terms give nonvanishing single cuts as well, we need to know them in more detail. 
We now compute the single cuts of some integrands in general form. The results will allow 
us to compute single cuts of spurious terms on one hand and the single cut of a total 
integrand expansion on the other. 



3. Single cut of general integrands 



In listing the single-cut results of general integrands, it is as convenient and more useful to 
list them in terms of general numerators and denominators. A general integrand will take 
the form 

nf=i(2fc-A) 



In,p{Ai, • • • , A p ; D , ■ ■ ■ , Ai-i) = T-r,,-! 



(3.1) 



since any appearance of the contraction k 2 in the numerator can be replaced by mf, in the 
single cut. Using ( [2 .7] ) in the limit u — > 0, we see that the single cut of Dq gives 



/•27T 

A Do [In, P ] =K 2 t p+1 / daG n , p (A u ...,A p -D ,...,D n ;a), 
Jo 



(3.2) 



where we have defined 
G ritP (A 1 , ...,A P ;D ,.. .,D n _i;a) = Ke l 



dz- 



nr=i(Ai^iA] 



The spinors (A| and |A] depend on z and z as follows, 

(X\ = (p\+z(ql \\] = \p]+z\q]. 



(3.3) 



(3.4) 



The single cut of is known once the integral G n ^ p has been computed. We now compute 
the single cut of I n ^ p for several values of (n,p). 
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Computation of Ad [Ii.i] 

It is easy to show that 



Gi,l(Ai; D ; a) = A 2 (p\A x \p] + A 3 e ia (q\A x \p] + ^ e -™ {p \A x \q] + ^-{q\A x \q}. (3.5) 



A 3 



A J 



Therefore, paying attention to the powers of t included in (3.2), we find the following values 
for the first two refinements of single-cut operators. 



A Do [I XjX (A x ;D )) =0. 

A Do [I X;X (A X ;D )) = ttK 2 A 4 (2A x • q)t 2 . 

Computation of Ad [I24] 

The A 2 t— part of the single cut is the A 2 /t part of G^i, which is 

A 2 Ai- q 
t K x -q 

The single cut of I2 1 is obtained after the trivial integration over a, giving 



(3.6) 



A Do [I 2>X (A X ;D ,D X )} 



Ai-g 
K x -q 



A 



D 



1 



(3.7) 



In the it — > limit, A/) [J2,i(Ai; Dq, D x )] is therefore proportional to the single cut of the 
tadpole. 

For later convenience we compute A^ [/2,i] and A m 2[/ 2 ,i] in the case where K 2 = 
and mg = m\ = m 2 . The logarithmic part of G% x is a lengthy expression. The logarithmic 
part of of the single cut of I 2 ,i is 



A Do [I 2A (A x ;D ,D x )} =2irK 2 t 



(p\K x \q]{q\A x \p] , (p\K x \p]{q\A x \q] 



+ 



+ 



(q\K x \q} 2 (q\Ki\q] 2 
A x \q](q\K x \p] 2(p\K x \q](q\A x \q](q\K x \p] <p|Ai|p] 



{q\K x \qf 
Y,C a (t)v%A llM . 



(q\KM 



(q\Ki\q] 



(3.8) 



The last line is simply an abbreviation for the expression, which will be convenient short- 
hand in one of the examples we give in Section [|. In that case where K 2 = and 
Trig = m? = m 2 , we also have 

A m 2 [I 2)X (A X ;D ,D X )} = A Do [I 2jX (A x ;D ,D x )} + A Dl [/ 2j i(A i; D , D x )\ = 0. (3.9) 

The single cut A/j 1 is computed redefining the loop momentum k — > k + K x , as described 
in Section |2[ 
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Computation of Ad [I2.2] 

The expression of G 2 ,2 is rather complicated and it is not shown here. We obtain Ad [12,2] 
(or Ad [-^2,2]) by expanding K 2 t 3 G2,2 in the limit A — > 00 and then integrating the 
A 2 t— terms (or A 4 t 2 — terms). The result is 



A Do [I 2 ^A 1 ,A 2 ;Do,D 1 )} = - h {A \' T f q) A Da 



Ai 



A Do [I 2 , 2 (A 1 ,A 2 ;D ,D 1 )\ = - 2vrif 2 A 4 {M ' q){M ' ' ; ' - 



(3.10) 
(3.11) 



For later convenience we compute A m 2/2,2 i n the case where K 2 = and = m 2 = m 2 , 

A m 2 [I 2t2 (A 1 ,A 2 ;D ,D 1 )] = A Do [I 2 , 2 (A U A 2 ; D , D x )] + A Dl [h^, A 2 ; D Q , D x )] 

= -{(2K 1 -A 1 )A Do [I 2A (A 2 ;D ,D 1 )] + (A 1 ^A 2 )}.(3.12) 

Computation of Ad [Is,i] 

The integrand K 2 t 2 G 3) i(Ai; Dq, D\, D2; a) does not contain A 2 t— terms. Therefore Ad [h,i] 
vanishes. 

A Aj [J 3j i(Ai;Z>o,-Di,^2)]=0. (3.13) 
This result confirms the absence of tadpole integrals in the reduction of I 3j \. 

Computation of Ad [13,2] 

The expression of G 3 2 (A X , A 2 ; Dq, D\,D 2 ; a) is rather complicated and is not shown. As in 
the previous cases, K 2 t 3 Gs !2 has to be expanded in the A — > 00 limit. We obtain Ad [Isp] 
by integrating the terms of 0(A 2 t), finding 



A A) [I 3>2 (A 1 ,A 2 ;D ,D 1 ,D 2 )} = [^[f^^ A A) 



(3.14) 



Computation of Ad [13,3] 

We obtain Ad [^3,3] by taking the A — > 00 limit of K 2 t 4 G 3t3 and integrating the A 2 t— terms, 
to get 

2 . (A, . n\( A n . n\( A„ . n\ I 1 

(3.15) 



A Do [I 32 (A 1 ,A 2 ,A 3 ;D Q ,D 1 ,D 2 )] = V f ^ 1 '^ 2 '^'^ A Dq 



4. Computation of tadpole coefficients 

In this section we show how the single cut allows the computation of tadpole coefficients 
in several examples of small integrands. 

We will compute a(0), the coefficient of the tadpole integral Aq^uiq), by the single cut 
operator A d defined in ( |2.23j ) . We look at one-loop integrands of the type 



n?=i(2fc-fr 



(4.1) 
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As described in the setup, when we cut the propagator Dq, we set Kq = 0. Any appearance 
of the contraction k 2 in the numerator should then immediately be replaced by jtiq, so that 
this form of the integrand is in fact general. 

For now, we assume that the Gram determinant of is nonvanishing and that the 
masses are non-degenerate. When this is not the case, further modifications are necessary, 
which we address in the following section. 

We will describe the derivation of the tadpole coefficients for the integrands with 
(n,p) = {(2,1), (2, 2), (3, 2), (3, 3)}. We have also verified the result for (n,p) = (4,3), but 
as this calculation does not involve any notable new features, we do not present it here. 
We have used FeynCalc [57] to check our results against those from Passarino-Veltman 
reduction. 

The idea underlying our procedure is to expand the integrand in such a way that 
spurious terms are easily recognized. The Ossola-Papadopoulos-Pittau (OPP) decomposi- 
tion [10] is well suited for this purpose. We take all the coefficients of both spurious and 
physical terms as unknowns. We then drop the physical terms except for the tadpole. Single 
cuts of all remaining terms are evaluated using the general results of the previous section. 
The tensor Y\i Ri is expanded in a basis constructed from fixed vectors; a convenient choice 
includes external momenta and orthogonal vectors, as in the OPP classification. Thanks 
to this expansion, the single-cut equation becomes a system of separate equations, which 
are the coefficients of independent tensors. In the following examples, we illustrate the 
derivation of the system and the solution of the tadpole coefficient. 

Tadpole coefficient of 12,1 

The starting point is the OPP decomposition of 12,1 given by 



'-2,1 



o(0) ~ 2k 

1 +Moi) 



D 



D D 1 



7 7 , s 2k 

+ 621 (01) 



D D 1 



~ , x 2k ■ n 



(4.2) 



Terms whose single cut contains no A 2 t contribution are included in "• • • " . These terms 
will be systematically neglected throughout this section. The momenta n, £7 and £% are 
defined [10] to satisfy the conditions 



Ki-n = Kf 



1 • £7 



Ki-£ 8 = 0, 



n 



It I* 



-K 



1 - 



£7 — I 



1 = 0. 



Applying the single cut operator A/j and using fl3.7|) , we get 



62i(01) -Q 4 )^+ (6o(01) -Q 2 



n' 



K x -q 



A 



Do 



D n 



(4.3) 



Here ai=i 1 ... ,4 are the coordinates of R\ in the basis {Ki,n, £7,^8}- Explicitly, they read as 
follows: 

Ri-K! Rx-n R 1 -£ 7 

"1 = 779 ) Q 2 = 772~ ' "3 = 7T2 - ' a 4 = 7T2 - • ( 4 - 4 J 



Kl ' 



Ki 
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From the relation (|4.3| ) we see our first example of a system of equations leading to 
the tadpole coefficient. Since q can be chosen arbitrarily, the expression inside the square 
brackets vanishes. This implies that each of the factors multiplying the basis vectors 
vanishes separately, giving four equations. 

In fact, we only need the first of these equations to get the tadpole coefficient; here we 
do not need to solve for any of the spurious coefficients. The result is 



a(0) + a.\ 







a(0) 



(4.5) 



This value is in accordance with the one obtained using the Passarino-Veltman decompo- 
sition. 

Tadpole coefficient of Z2 2 

For 2^,2, again we use the OPP expansion, keeping only the terms with a non- vanishing 
single cut of D$. These are the tadpole along with spurious terms, 

o(0) r , ,2k-£ 7 ~ , 2k ■ £ 8 ~ 2k ■ n 
^2,2 = 6ll(01) 7^ + 62i(01) — -2 + 6o(01) 



^0 
+ 6i2(01) 



D D 1 
(2k ■ ^ 2 



D D 1 

\2 



D D 1 



;r + 622(01) + moi) 



D D 1 



D D 1 
(2k ■ l 7 )(2k ■ n) 
D0D1 



~ , . (2k ■ £ 8 )(2k ■ n) ~ , , 
+ 602(01) ^ ^ " + ft oo(01) 



D Dt 
(2k-K 1 f -Ak 2 Kl 
3A)#i 



(2k ■ nf 
D D l 



+ 



(4.6) 



We compute the single cut of both sides of (|4.6|) using (|3.7|) and ( 3. 10| ) - The outcome is 
= 



O(0) + ai + 6oo^l) \ K » K „ + _ \ ^ + /^ 2(01) _ \ ^ 



+ (6 22 (01) - 04) iffi + 

'621(01) 



2a 7 Kffi 



+ 



fx 



2a 8 K»l\ + 6oi(01) - 2ai 3 n»t 7 



+ 602(01) -2a u 



A 



D 



D 



(K x -qf 



fx- 



(4.7) 



Here f\ is defined according to equation ( |2.12| ), while Q!i=i ) ... ) i6 are obtained from the 
following decomposition of R\ in a basis of independent tensors: 



R%B% = axK^Kl + a 2 n^n v + ^ [04-4^] + a 5 (-^fn^ + i^n") + a 6 (Kfn' y - 

8 

+ Y, + + «*+2(^r - ^)] + "11 (4^8 - W + «125m- 



i=7 
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(4.8) 



i=7 



The condition ( |4.7[ ) is fullfilled for any lightlike g only if the second rank tensor inside the 
square brackets vanishes. The first independent tensor includes the tadpole coefficient, but 
also the spurious coefficient 600 (01). So we need the second independent tensor as well, but 
no others. In particular, 



f a(0) 6qo(01) 
k 6oo(01) - a 2 = 







a(0) = -/i ai + 



Q'2 



(4.9) 



Using the explicit expression of a\ and «2; 

3(#i • R^Kt -R 2 ) (i?i-i? 2 ) (n-/2i)(n-i?2) 



a 2 



2(^x 2 ) 2 



2^ 2 



2(1C? 



2\2 



2(^ 2 ) 2 

we get the value of the tadpole coefficient 

/1 



2^ 2 



+ 



2\2 



a(0) 



^ {R 1 -R 2 )-A{K X -R 1 ){K 1 -R 2 ) . 



(4.10) 



3(^ 2 ) 2 

Tadpole coefficient of Z32 

We now find it convenient to vary the OPP decomposition slightly. Our single-cut decom- 
position of 2^2 reads as follows, 



!-3,2 



q(Q) 
d 

2 



+ ^c i2 (012) 



(2Jfc • £ 



2 2 



i+2j 



1=1 



DqD\D 2 



E^o(Oi 



i=l 



2/c • if 3 _ 
DoDi 



3=1 i=i 
• K 2 2A; • if, 
if 2 DoA 



2^+2 



+ 



As in OPP, ^3 and £4 are lightlike momenta such that 



(4.11) 



(4.12) 



q = Ki-£ } = 0, Vie {1,2}, je{3,4}. 

The evaluation of the single cut is performed using equations Q, (|37To|) and The 
result can be written as follows, 







- (6o(01) + 7 2 2 ) K>KZ + (a(0) - 27+ + 6(01)^p + 6(02)^^) J£J 



^ (5,2(012) - 7W2) ^+2^+2 - 6n(01) + 2 72 + 3 K!& - 6 n (02) + 2 7 + ) 



i=i 



621(01) + 2 72 + 4 ) K»tl - (621(02) + 2 7 + ) - (6o(02) + 7^ ) 



x A 



Do 



Da 



(4.13) 
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The coefficients 7^' are obtained from the decomposition of RiR 2 in a basis of independent 
tensors, 

2 

i=l 



+ E E [-rSW? + «r?) + 7«(^ - *r*?) 



i=l j=3 



+ 7 3 4 



3^4 



(4.14) 



Eq ( 4.13 ) implies that the following conditions have to be fulfilled, 

6o(01) +7 2 ° 2 = 0, 

4- ~ , s K 1 ■ K 2 r , , Ki ■ K 2 
a(0) -2 7 ++ bo(01)^2^ + 6 (02)- 



6 (02) + 7 °! = 0. 



Kl 



0. 



(4.15) 



The system ( [4.15 ) together with the explicit expressions of 7^, 722 and 7^, gives the 
tadpole coefficient of Z32, 



a(U) - 2 7l2 + 722— + 7u-^- " K 2 K 2(( Kl . K2 )2_ K 2 K 2y 



(4.16) 



The factors bij appearing in ( 4.16 ) are defined as follows, 
612 = 6 2 i = -KlKl 6n = K\(K X ■ K 2 ), 



>22 



K\(K X -K 2 ) 



The decomposition (4.11) relies on the particular structure of the numerator of 2^2- 
Being more general, the OPP decomposition does not take advantage of the knowledge of 
the numerator of 2^2- As a consequence, new spurious terms enter. They are of the type 



(2k-P)(2k-Q) 



2k ■ P 



(4.17) 



with P,Q G {Ki, K 2 ,£%, £4} and i = 1,2. The coefficients of the terms ( 4.17| ) vanish. 
This can be explicitly shown using the operator Ae> defined in (2.24), which selects the 
A 4 i 4 -enhanced terms of the single cut of 2^2. 

Tadpole coefficient of Z33 

This is the last example of a tadpole coefficient that we will describe. We have done the 
analogous calculation for 14,3, but it does not introduce any notable new features. 

Here again, although we could have started with the full OPP expansion or other vari- 
ations, we can simplify the calculation by the particular expansion of 2^3 in the following 
tensor integrands: 



-3,3 



«(0) 
Do 



+E* 

i=3 



(2k 



+ EE c » 

i=3 j=\ 



D DiD 2 ' 



(2k ■ £i)(2k ■ £i)(2k ■ Kj) 
D DiD 2 



- 13 - 



4 2 2 

+ EEE 

Hjk 

i=3 j=l k=j 



(2k ■ £i)(2k ■ Kj)(2k ■ K k ) 



D D 1 D 2 



E Ciij I fj 

ij=i \ 



(2k ■ Kj) 



+ 



sp. 



(2k • Kj) 2 



(4.18) 



sp. 



Here we are using the definition of £3 and £4 given in ( |4.12 ). The operator [• • -] sp . selects 
the spurious parts of its argument. The explicit expressions of the spurious part can be 



read from equations ( f4.2|) and ( [4.6|) . The single cut of both sides of equation ( 4.18| ) is 
computed using ( |3.7| ), ( p. 10 ), and ( 3.15 ). Using the decomposition of R^R 2 R 3 in a basis 
of independent tensors, 



R»R v 2 Rl= E ^K^KZ + ^^^ + ^^o^sT' + d^W 

i,j,k=l i=3 i=3 

2 

+ <Wf ~ m) ) + E { a mK?3 ua + c&Krgr + 

i=l 

2 

+ E (Voo*? ^ - m%) + «o,o^r - mi) + «o ^r - «) 

i=l 
2 4 

+ E E {^k K i K ^l + + a kij KVK?t^) 

(4.19) 

and equating the single cuts on both sides of equation fl4.18j ), we get the following relation: 

2 

»=E 

m=l 
2 

+ E m(c) K»K\Kl + 53 - m m; m(5) K^K^K^ 
i=l 

+ ■^-K^Kl^ ( o(0) - E 



j J=l fc=3 
2 4 

+ E E + m K i*ffi + > 

i=l j=3 



E e ^ W + E E E W + E E 

i=3 i=3 j=l fe=_; 



i=3 fe=l 



Hjk 



^fm 9m 3— m; m(c) 



A 



Do 



(K x -q)(K 2 - q )(K m -q) 



(4.20) 



The coefficients in the above relation have been defined for convenience of displaying the 
independent tensors. We use the following abbreviations (i ^ j ^ k ^ i), 



C-iii — Cm OLiU 



Ciik — Cjjfc (ctkii ~i~ OLiki ~i~ Qiiife) 

Cijk = Cijk ~ [ctkji + otkij + ctikj + (1 - 5jk)(otjki + ctjik + a ijk )} , (4.21) 
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and the following definition of the totally symmetric coefficient dijk, 



fi 



Ki ■ K 3 „i f 3 _i ( K 2 A (K 3 -i ■ Ki 



K 2 
A 3-i 



K 2 



fi + f: 



The gij. m (c) are given by 



1 



^1^2 



(if \K\ - \{K X ■ K 2 f) cjjj + 



K x -K 2 



3(K 2 ) 2 

Since q is arbitrary, we are led to the following relations. 

' 311; 2(c) = 
322; l(c) = 

/l3ll; l(c) + /2 312; 2(c) = 
/l321; l(c) + h 322; 2(c) = 



2 



(4.22) 



(4.23) 



(4.24) 



a (°) ~ ^2 a ijk d ijk + fl 312; l(c) + /2 321; 2(c) = 

which uniquely fix a(0), 

2 

a(0) = a ijk dij k . (4.25) 

2 J,fc=l 

We have observed that a different choice of the original integrand expansion, rather 
than ( 4.18Q , can lead to a larger linear system which is not completely solvable, yet there 
is still a unique solution for a(0). The underlying phenomenon is that different spurious 
terms can have the same single cut, so these terms should really be grouped together in 
the expansion. 



5. Massless external legs 

As in other techniques such as Passarino-Veltman [21-25, 29, 30, 35, 58, 59] or Ossola- 
Papadopoulos-Pittau [10,50], we must modify our algorithm in the case of vanishing Gram 
determinants. The most immediate case is the presence of a massless external leg in tensor 
bubble integrals. In this section we will focus on this class of integrals in the case where 
the internal masses are the same, uiq = m 2 = m 2 . 

New features appear in the computation of a tot , the coefficient of Ao(m 2 ). First of all, 
if K\ is light-like, it is not possible to complete a basis defining n,£j,£s. A suitable basis 
is given by {K\, K 2 , £7, is}, which is composed of four light-like momenta such that 

KfK 2 ^0^£ 7 -£ 8 , K % -1 3 =^ ViG{l,2},jG{7,8}. (5.1) 
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Secondly, the coefficients of 1/Dq and 1/D\ contribute to a tot since 

Moreover the scalar bubble and the tadpole are connected, 

A {m 2 ) 



B (0,m 2 ,m 2 ^ 



1. 



(5.2) 



(5.3) 



so the bubble coefficient contributes to the total tadpole coefficient. Tensor integrals con- 
tracted with K 2 are no longer spurious terms [50] and they contribute to the total tadpole 
coefficient. This can be easily understood looking at their explicit expression [35], 



/ 



H^ = i M 0,™W)(2AV* 1 )- I ''- 4 ° (m2) 
(2k • K 2 ) 2 



<Fk 
d 4 k 



B u (0,m z ,m z )^K 2 -K h 
L>oL>i 6 



2 1 {Mm 2 ' 



lj (2K 2 -K X ), 

- (2K 2 ■ K x ) 2 . (5.4) 



The tadpole coefficient a tot is obtained by summing the aforementioned contributions. 
These contributions can be obtained by applying the single cut operators defined in ( 2.23; ) , 
( |2~24|) , and ( gjg ). 



Tadpole coefficient of 12,1 
The full integrand is expanded as 

a(0) a(l) b 
-L-2,1 = — r — h ^ ^ + • 



D x D D 1 



Sd(oi)^^+][>(oi) 

L> L>1 



i=l 



{2k ■ e i+6 ) 

D D 1 



+ 



The terms denoted by "• • • " are those such that A m 2 [• 
tadpole coefficient is given by 



o(0) + o(l) + \ + • K X ). 



(5.5) 
0. The 

(5.6) 



The sum a(0) + a(l) can be obtained by applying A m 2. Using equation ( |3. 9[ ) we see that 
only the tadpole terms survive, so 



a(0) + a(l) = 0. 



(5.7) 



The coefficient 6o(01) is obtained by cutting the propagator Dq and selecting the A 2 t terms. 
The outcome is 



(o(0) + ai)K? + (S (01) - a 2 )K% 



+ ^(M01)-a i+2 W +6 



Ad 



(5- 
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The parameters ai ... 4 are the coordinates of R\ in the basis {Kx,K2,l7,ls}, which read 
as follows: 



R x -K 2 Rx ■ K x 

Q l = "77 TTi a 2 



R 



Rl 



Since q is arbitrary, each coefficient of the momenta appearing in ( |5.8D has to vanish. In 
particular, 



6o(01) = a 2 



(5.10) 



and 6ji(01) = ai +2 - Finally, the bubble coefficient b can be obtained by using the operator 
Ad defined in (2.24). Using equation (3.8), we get 

r r 2 

Q=J2 dt c a(t)v a ^ - a x K^ + (6o(01) - a 2 )K% + ^ (&ii(01) - a l+2 ) 



l i+6 



+ 6 A 



D £>i. 



D Dx 



b A 



Do 



1 



0. 



The tadpole coemcient is obtained from equation (| 
and Q5.ll ), and it is given by 

(Ri ■ Kx) 



(5.11) 

using the results Q, (fig ), 

(5.12) 



as expected from explicit reduction. 



Tadpole coefficient of 

Our final example is given to display the flexibility of single cut operations and the degrees 
of information available from cutting the same propagator. We use several different refine- 
ments of the single cut, picking out terms with different dependence on A and t, in order 
to collect the subset of the information required for the tadpole coefficient. 
The integrand 12,2 is decomposed as follows, 

2 



'-2,2 



a(0) o(l) 



D Dx 

+ £ (moi) 

1=1 ^ 



i=i 



2/fc • iT, 



+ a».+2(0) 



2Jfe •£ 



i+6 



(2k-K 2 )(2k-l i+ _^ + ~ (2k 



+ 6oo(01) 1 n n V +••• • 

The terms denoted by "• • • " are not explicitly shown since A r 
0. The tadpole coefficient of Aq(ui 2 ) is given by 



(5.13) 



m 



b_ M 6oo(01) , 
2 



Otot = o(0) + o(l) + £2 + ^g^(2K 2 • i^) 2 . 



(5.14) 
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In the computation we will take advantage of the following decomposition of RiR 2 m 
independent tensors, 

2 8 



8=1 



i=7 



8=7 

8 

+ a 12 g^ + J] [oi+eC^ + K^) + a i+8 (i^ - K%t%)\ . (5.15) 



i=7 



We apply the single cut and we select the A 2 t terms using the operator A m 2 and equation 
'a(O) +o(l) -4ai2 



dUD. 



2K X ■ K 2 
+ ^(6oi(01) - 2a m2 )C 



2a 5 + 2(6 o(01)-a2)^ 2 At 



+6 



i=l 



2(iTi • K 2 )g t 



■A 



Do 



1 



'oj 



Equation ( [5. 16 ) fixes the values of a(0) + a(l) and of 6oo(01) to be 

o(0) + a(l) = 4ai2 + 2a 5 (2^i • K 2 ), 



(5.16) 



(5.17) 
(5.18) 



^00 (01) = «2- 

The value of 6oj(01) is fixed as well, to be 

6w(01) = 2a i+12 . (5.19) 

The value of 6^(01) is read off from the A 4 t 4 -enhanced terms of the single cut of X 2;2 . 
These terms are selected using the operator Ajj and using equations (^) and ( 3. . 
The outcome can be written as follows, 

i 







- (5i(0) + a x )K^ - (a 2 (0) + 2a 5 )K?K% - ^(a,(0) + 2^+4)^+4 

i=3 

2 2 

+ E(M01) - 2a i+12 )K%£» +6 + J2(b i2 (01) - a l+2 )t l+& i 
-i=i 



'V 

i+6 



+ (6oo(01) - a 2 )K%K> 
6i2(01) = a i+2 . 



2q li q v 
Kt-q 



ttK 2 AH 2 



(5.20) 

The logarithmically enhanced part of the single cut allows the computation of the bubble 
coefficient b. This is achieved by applying the Ad operator defined in ( [2.24 ). Using 
equations ( 5.19 ) and ( 5.20| ), we find 



= b A 



Do 



D D 1 



4a 12 m 2 A Do 



b = 4m 2 a 12 . (5.21) 
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Here we have used the property Ad I2,2(Ki, A 2 ; Dq,D{) = 0, which holds in the case f\ = 
0. The tadpole coefficient in Equation ( |5,14| ) is then obtained using equations ( 5.17 ), ( 5.1§| ) 
and C|5T2lD , 

a tot = 2 (4ai2 + a 5 (2Ai • K 2 )) + a 2 ^—^ = 2R 1 -R 2 -\ ^—^ . (5.22) 

In obtaining the result (|5.22| ), we used some explicit values of the coefficients in the expan- 
sion ( 5.15| ), 

4a 12 + a 5 (2^i • K 2 ) = (R x • R 2 ), a 2 = (5.23) 

(K 1 ■ K 2 y 

6. Conclusions and discussion 

We have seen how single cut integrals can distinguish scalar boxes, triangles, bubbles and 
especially tadpoles. We have outlined a strategy to find the tadpole coefficients applying 
single cuts. In particular we have used the fact that the null vector q was chosen arbitrarily 
in order to establish the independence of contractions with it and derive a system of linear 
equations. The tadpole coefficient may been obtained as a solution of this system. 

The general application of our procedure requires further investigation. Since single 
cuts of spurious terms do not vanish, we need to understand the single cuts of general 
integrands in order to isolate tadpole coefficients. We have made a start in this direction in 
Section ||. It would be very interesting to generalize this analysis. The simple expressions 
of the terms proportional to A 2 t lead us to speculate that further results will be similarly 
simple. Perhaps there is even a more direct way to derive these terms in particular. 

Although in the examples in Section [|, the system of equations derived from single 
cuts was sufficient to uniquely determine the tadpole coefficient, we need more general 
information about the system to be sure that it will always work. The rank of the matrix 
associated to the system of linear equation could be insufficient to find a unique solution 
for the tadpole coefficient. A possible concern is that q is null and thus has fewer degrees 
of freedom than a generic vector. Here again, study of single cuts of general integrands 
could illuminate the role played by q and perhaps resolve our concern about the general 
solvability of the linear system. 

In Passarino-Veltman reduction, a single system of equations is solved for all coefficients 
of master integrals, and no spurious terms are involved. It is also unnecessary to reduce all 
the way to the tadpole level. In our procedure, by contrast, we actually solve two systems 
of equations: the first one is the expansion of the tensor numerator in a suitable basis of 
independent tensors. However, the second system leading to the tadpole coefficient is then 
quite simple and usually reducible. We have assumed that the coefficients of other master 
integrals would be obtained by (generalized) unitarity methods, so the total calculation 
is in fact quite long. We cannot claim that our method for getting tadpole coefficients 
will be efficient. Nevertheless, we have found it illuminating to investigate the single cut 
integral formalism in general. It would be very interesting to further probe the analytic 
structure of amplitudes, for example by developing a D-dimensional extension for studying 
their rational parts along the lines of [17]. 
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A. Phase space parametrization 

In this appendix we describe some details of the evaluation of the single cut. The starting 
point is the integral 

(PhS+dtft g, (A.l) 



where g is a general integrand. Given an arbitrary four-vector K, such that K 2 ^ and 
its energy component is positive, Kq > 0, we construct a pair of light-like momenta p and 
q such that p + q = K, and two more momenta ei, €2 defined as 

ei = ~(<g|Y*|p] + <Pl7"l?]), f 2 = ^((g| 7 M N - <P|7"M 
The loop momentum t\ can be decomposed in the basis {p, q, ei, 62}, 

h = tip + aq + e\x - e 2 y), 



and the integral ( A.l ) can be expressed in terms of the coordinates (t, a, x, y). The Jacobian 
of this reparametrization reads as follows: 

<T _ I /3_ r ^n v e p f : (T \ - 2 ( R2 } 2 1+31 

J — I 1 E y,vpoV Q e i e 2 I — ^ r I ' 
where in the last step the relation 

U e^a^Cpdv = -2{a\bcd\a] + (2a • 6) (2c • d) - (2a • c)(2b ■ d) + (2a • d)(2b ■ c), 
has been used [60]. The integral ( |A.l| ) becomes 



d 4 e 1 5 + (£ 2 1 ) g = I dtdxdydaJ5 + (t 2 (a - x 2 - y 2 )) g 

00 f K 2 t 

dt / dx dy g. (A. 2) 

Jr2 2 

This integral is computed using the generalized Cauchy formula, 

2iri V Res{F(z,z),Zj} = <k F(z, z)dz - [ ^dz A dz. (A.3) 
Jt?,. Jds Js oz 



poles Zj 



Considering g as a function of (z, z) = (x + iy, x — iy), we construct its primitive G(z, z) 
with respect to z. Choosing S to be the complex plane C, we get 

/ dx dy g = — / -—dzf\dz = — / G(z,z)dz — tt > Kes{G(z,z),Zj}. (A. 4) 
Jr2 2z J c dz 2i J dc j-f 



poles Zj 
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In practice we regularize the divergences of the integral ( |A,4j ) taking instead of C a disk 
D in the complex plane enclosing all poles of G(z, z) viewed as a function of z. In polar 
coordinates the disk is parametrized as follows, 



z = re 



D = {(r, a) | < r < A, < a < 2tt}. 



Under these assumptions, equation (A^4) becomes 



1 f 2?r 

dxdyg= lim - / Ae ia G (Ae ia , Ae~ ia ) da - vr V Res{G(z, z), zA. (A.5) 

A^OO 2 In ' Z — ' 

u poles Zj 



In equations (A.2) and ( |A,4| ) the correct prefactors are kept, even though in unitarity 
methods these constant factors cancel out and can be neglected. 2 
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